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We describe a many-body quantum system which can be made to quantum compute by the 
adiabatic application of a large applied field to the system. Prior to the application of the field 
quantum information is localized on one boundary of the device, and after the application of the 
field this information has propagated to the other side of the device with a quantum circuit applied 
to the information. The applied circuit depends on the many-body Hamiltonian of the material, 
and the computation takes place in a degenerate ground space with symmetry-protected topological 
order. Such adiabatic quantum transistors are universal adiabatic quantum computing devices which 
have the added benefit of being modular. Here we describe this model, provide arguments for why it 
is an efficient model of quantum computing, and examine these many-body systems in the presence 
of a noisy environment. 



I. INTRODUCTION 

The invention of the transistor [1] was a watershed 
moment in the history of computing: it provided a logic 
element that was naturally robust to noise and error. 
Quantum computers offer the potential to exponentially 
speed up some computational problems (notably factor- 
ing [2]), but have not been built in large part because 
quantum information is notoriously fragile and quickly 
becomes classical information in the presence of noise. In 
theory, the quantum threshold theorem [3-5] asserts that 
these difficulties can be circumvented, but in practice the 
requirements of this theorem are daunting. Here we out- 
line a novel method for building a fault-tolerant quantum 
computer that much more closely mimics the classical 
transistor. In particular, we show how a suitably engi- 
neered material can quantum compute by the simple ap- 
plication of an external field to the sample. Applying the 
field adiabatically causes quantum information to spa- 
tially propagate across the device at the same time that 
a quantum computation (quantum circuit) is enacted on 
the quantum information, and this in turn allows us to 
design clocked quantum computing architectures similar 
in control requirements to modern classical computers. 
While we will not be able to rigorously show that our 
adiabatic quantum transistors are fault-tolerant devices 
like classical transistors, we present analytical and nu- 
merical arguments as to why these transistors could be 
tolerant to errors and thus a true quantum analog of the 
classical transistor. 

The standard operating model for a quantum com- 
puter is called the quantum circuit model [6]. In this 
model, one begins with a system of initialized two-level 
quantum systems (qubits), applies a temporal sequence 
of one and two-qubit gates enacting a circuit, and fi- 
nally performs a measurement (readout) of the qubits. 
In contrast with modern classical computers where in- 
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formation propagates spatially — advancing roughly one 
step across the computer chip at each rise and fall of a 
clock voltage — most proposed implementations of quan- 
tum computers fix the information spatially and bring 
the computational operations to the data. Many re- 
searchers have noted this difference, and a variety of 
quantum computing (QC) models were subsequently de- 
veloped where quantum information propagates spatially 
during a quantum computation; examples include spin- 
wave models where quantum information moves ballis- 
tically down a quantum wire [7], linear-optics QC [8], 
one-way QC where simple sequential measurements push 
the quantum information across the system [9] , and some 
universal adiabatic QC models which create a superposi- 
tion of computational states spread across the device [10- 
12]. However none of these constructions yield architec- 
tures that mimic modern synchronous sequential com- 
puter chips. In these chips the rise or fall of a global 
voltage is the trigger that causes the information in the 
device to move spatially across the device. While the 
control requirements for the clock signal in synchronous 
sequential logic devices are by no means trivial, they do 
not require the precise control, measurement, and tim- 
ing that makes quantum computers notoriously difficult 
to build. Here we introduce a new way to mimic these 
classical control requirements that allows us to propose 
synchronous sequential fault-tolerant QC architectures. 
The key to our construction is a novel type of quantum 
gadget that mimics the role classical transistors play in 
modern computers, and that we hence label an adiabatic 
quantum transistor. While the outline we give for an 
adiabatic quantum transistor is still very far from ex- 
perimental realization, the novel quantum computational 
matter that such a device represents opens a new and po- 
tentially promising path for the construction of a large 
scale quantum computer. 

An outline of the paper is as follows. In Sec. II we in- 
troduce a model of an adiabatic quantum transistor. The 
speed at which we can operate a quantum transistor is re- 
lated to the minimal energy gap between the ground state 
manifold and first excited state during the application of 
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the applied field. In Sec. II A we rigorously prove that 
if the adiabatic quantum transistor enacts a single-qubit 
quantum circuit with only identity gates, then the energy 
gap is such that the adiabatic quantum transistor can be 
enacted efficiently. In Sec. II B we provide strong nu- 
merical evidence from matrix product state simulations 
that the same efficiency holds when the adiabatic quan- 
tum transistor enacts an arbitrary single-qubit quantum 
circuit. Next in Sec. II C we provide (weaker) numerical 
evidence that the same results holds for our adiabatic 
quantum transistors for quantum circuits involving more 
than one qubit. Having given evidence that the adiabatic 
quantum transistor model is an efficient model of quan- 
tum computation in an ideal word in which the system 
is isolated from its environment, in Sec. Ill we turn to 
the question of whether the model can be made fault- 
tolerant. We give arguments that if quantum transistors 
are configured to execute fault-tolerant quantum circuits 
then the tolerance of these circuits to errors is conveyed 
to our model. While we cannot show that our model has 
a threshold theorem associated with it, we can give phys- 
ical reasons for why the model will have a threshold. In 
Sec. IV we discuss how the unrealistic four-body Hamil- 
tonian we use in our constructions can be implemented 
using a Hamiltonian with only two-body interactions, via 
perturbation theory gadgets, and discuss the effect that 
these gadgets have on the arguments in the previous sec- 
tions. In Sec. V we present a variety of constructions for 
quantum adiabatic transistors that are likely to be use- 
ful, including systems that can be used to spatially route 
quantum information and to perform measurements. Fi- 
nally in Sec. VI we conclude and list important open 
questions for the adiabatic quantum transistor model. 



II. THE ADIABATIC QUANTUM TRANSISTOR 
MODEL 



Our proposed adiabatic quantum transistor is a device 
that operates in a manner similar to a classical logical 
element such as a MOSFET. See Fig. 1 for the analogy 
with a classical transistor and Fig. 2 for a diagram of 
how an adiabatic quantum transistor would work under 
this analogy. Like a classical transistor this device is 
made to quantum compute via the application of an ap- 
plied field to the device. Prior to the application of the 
field, the quantum information is localized to one side of 
the device, and the system is in its ground state. Af- 
ter the application of the field, the quantum information 
is on the opposite side of the device, but with a quan- 
tum circuit applied to this information, and the system 
remains in its ground state. The exact circuit applied 
depends on the microscopic details of the system — we 
propose that each such quantum transistor be made to 
implement a fault-tolerant quantum gate or a gate that 
can be used to classically steer the quantum information 
across the device. During the application of the field, 
the system remains in its ground state, however the en- 
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FIG. 1. Graphical depiction of the operation of a classical 
transistor. In a classical transistor such as the n-channel 
MOSFET shown, the application of a voltage at the gate 
input transforms the semiconductor in a manner such that 
a channel opens between drain and source. Note that this 
transition is switched by the application of an electric field. 
Traditionally we think of this action as a switch, but there 
is another possible interpretation. In particular, the appli- 
cation of the applied electric field has the effect of copying 
the voltage at the source to the drain, which we can inter- 
pret as an identity logical gate (a logical gate that takes in 
a single bit and outputs this bit unchanged) on the classical 
information encoded as a voltage. Quantum information can- 
not be copied [14], so if we are to image a suitable analogy of 
the transistor in the quantum world, the quantum informa- 
tion must move from the source to drain. This is what our 
proposed adiabatic quantum transistors achieve. 



ergy gap to the first excited state gets smaller. Thus 
the device functions as an instantiation of a universal 
adiabatic quantum computer [10-12] with, however, the 
important differences that the ground state of the sys- 
tem can be degenerate (and thus the model could also be 
considered an open loop holonomic computation [13]), 
and, more importantly, that the information propagates 
spatially (thus the transistors are modular.) This allows 
us to design clocked architectures that very closely re- 
semble today's classical synchronous digital computers 
(see Fig. 3). Here we introduce the details of adiabatic 
quantum transistors and present theoretical arguments 
and evidence from simulations that the energy gap in the 
system shrinks sufficiently slowly so as to allow adiabatic 
quantum computing. We further present arguments for 
how our construction can be made fault-tolerant. This 
results in an implementation of quantum computers in 
which precision gates, preparations, and measurements 
are replaced by sufficiently slow and smooth application 
of fields, and the degree to which one can suitably engi- 
neer a many-body quantum system. 

Our starting point is the recent work showing how it 
is possible to use piecewise adiabatic deformation of a 
many-body interacting quantum system [15-18] to per- 
form a quantum computation. In [16] we considered a 
certain many-body interacting quantum system in its 
ground state whose Hamiltonian we adiabatically de- 
formed by turning on a strong field at the border of the 
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FIG. 2. Sketch of how an adiabatic quantum transistor be- 
haves. With no applied field, the quantum information is 
on one side of the device (shaded grey, left) and after the 
adiabatic application of the field the quantum information 
has propagated to the other side of the device (shaded grey, 
right) with a unitary U applied to the original information. 
The strength of the applied field is initially zero and then is 
ramped up to a non-zero value, during which time the quan- 
tum information is not localized to either end of the device. 
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FIG. 3. Using quantum transistors one can design a syn- 
chronous architecture as shown. A router takes classical infor- 
mation about what circuit should be executed (the quantum 
program) and routes it to the appropriate action unit. These 
units are made up of a collection of quantum transistors that 
can be used to implement their given functions. This infor- 
mation is then directed to the router where another step in 
the circuit begins by reading the next program action. 



device and then slowly propagating this field across the 
material; this propagates quantum information through 
the device in a piecewise fashion, riding the front of the 
applied field. A major drawback of this method is that 
it requires that one precisely control a microscopic ap- 
plied field which is turned on in a step-by-step process. 
This motivates the question: if we take the construction 
in [16] and, rather than turning on the field in a piece- 



wise fashion, we instead turn on the field simultaneously 
across the entire device, then does the device still per- 
form the desired quantum computation? Here we answer 
this question affirmatively and show how it leads to our 
adiabatic quantum transistor. 

The many-body quantum system with the remarkable 
properties described in [16] is described by a twisted ver- 
sion of the Hamiltonian associated with cluster states. 
Cluster states are the entangled states originally used 
to perform measurement-based quantum computation [9, 
19]. To define the Hamiltonian of our system, we need a 
graph (with vertices V and edge set E) where each ver- 
tex, v, is labeled by an angle V and is associated with 
a qubit. We define the twisted cluster- state Hamiltonian 
to be 
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where [P] v denotes the operator P acting on the qubit 
at vertex v, X, Y, and Z are Pauli operators, and we 
choose units so that the coupling constant is unity. The 
ground state of this Hamiltonian is a cluster state in a 
locally rotated frame [16]. Note that the above Hamilto- 
nian contains interactions that are physically unrealistic 
because they involve more than two qubit s; fortunately, 
it turns out that we can use perturbation theory gadgets 
to obtain an effective Hamiltonian of this form inside 
the low-energy sector of a Hamiltonian containing only 
two-qubit interactions [15, 20, 21] (see Sec. IV). We will 
assume here the Hamiltonian in Eq. (1) and return to 
implementation in two-qubit interactions later. 

The results in Ref. [16] show that one can take a quan- 
tum circuit and, using the recipe shown in Fig. 4, map 
it to a twisted cluster-state Hamiltonian with the prop- 
erty that adiabatically turning on fields along the x di- 
rection (— X v on vertex v) and turning off relevant non- 
commuting terms in He pushes information through the 
device in such a way as to implement the computation 
defined by the quantum circuit. In the analysis of [16] 
one can separate out the rigorous argument that the com- 
putation is performed from whether this procedure can 
be performed adiabatically. Thus, if we consider the case 
of turning on all of the fields simultaneously rather than 
piecewise, we see that our analysis of the fact that a quan- 
tum computation is performed carries directly over from 
[16]. By contrast, although the energy gap in the com- 
putation is independent of the size of the computation in 
the piecewise construction of [16], when one turns on the 
field at all locations simultaneously this will no longer 
be true. Therefore, in order to show that turning on the 
field over the entire system still causes our device to ex- 
ecute the desired computation, we must show that the 
inverse of the energy gap between the ground state and 
the first excited states grows at most polynomially in the 
size of the quantum circuit being enacted. If indeed the 
inverse gap grows this slowly then the adiabatic theorem 
(for example [22]) guarantees that the system will remain 
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FIG. 4. Dictionary for translating a quantum circuit into 
the cluster state Hamiltonian (see one-way QC [9].) First, 
express the quantum circuit in the universal gate set above 
where Rq — exp(— iOZ/2), if is a Hadamard, and the two- 
qubit gate is a controlled-phase gate, along with preparations 
of the +1 eigenvalue of A, |+). Then convert this circuit into 
a graph and corresponding twisted cluster-state Hamiltonian 
as diagramed in the dictionary, merging white output qubits 
to shaded input qubits in the circuit. 



in its ground state throughout the evolution conditional 
on turning on the field over a time scale polynomial in 
the size of the quantum circuit. 

To summarize, quantum circuits correspond to via 
Fig. 4 to twisted cluster state Hamiltonians, and the 
minimal spectral gap of such a Hamiltonian determines 
how long it takes, via adiabatic deformation, to achieve 
a high-fidelity implementation of the circuit. 



A. A One Dimensional Wire 

We begin by analyzing the simplest instantiation of our 
model: n qubits on a line with no twists in the Hamil- 
tonian, 6 V = for all v G V. In this case we evolve the 
system according to the time varying evolution 

n-2 



H(t) = - f(s)[y j [Z] l [X} l+1 [Z} l+2 - [Z] n _i[X] n 

n-l 

-sWE(4 (2) 



where / and g are suitably smoothly varying envelopes 
which satisfy /(0) = g(l) = 1 and /(l) = g(0) = 1, and 
s — |? is a scaled time. For simplicity, in this subsection, 
we will assume f(s) = 1 — s — l—g(s). Note that we have 
also turned off the cluster state term in Eq. (1) corre- 
sponding to the first qubit. This implies that the ground 
state of the Hamiltonian is two-fold degenerate and thus 
encodes a qubit [16]. Initially this qubit is localized on 



the first two qubits of the chain. At the end of the adia- 
batic evolution the qubit will be at the end of the chain 
on the last qubit, with a Hadamard gate applied to the 
qubit if the chain has even length. We can now ask the 
question, what is the minimum energy gap for the above 
evolution, and how does it scale with the length of the 
chain? If the energy gap scales as an inverse polynomial 
in the length of the chain, then the adiabatic theorem 
tells us that propagating the qubit down this length n 
chain can be done efficiently (i.e., in a time polynomial 
in n). 

To answer this question, we use the equivalence of this 
model to two uncoupled transverse Ising models [23] (see 
also [24] ) . These systems can be exactly diagonalized by 
a transformation of the spin model to a model with non- 
interacting fermions. Following the work of Doherty and 
Bartlett[23], but with additional attention to the bound- 
ary terms, suppose that one defines a code with the fol- 
lowing (n — 1) encoded Pauli operators 



[TO- 

[x] 2 [xu- 



[X]j if j is odd, 
[X]j if j is even . 



and 



[Z]j[X) j+1 [Z) j+2 if j<n-l, 
[Z]„_i[*]„ ifj = n-l. 



(3) 



(4) 



Then, given these encoded operators, we can then express 
the Hamiltonian in Eq. (2) as the union of two uncoupled 
transverse Ising models acting separately on the even and 
odd qubits, with an additional boundary term: 

n—l / n — 3 \ 

H(s) = -(1 - s) Zi ~ s ( X x + X 2 + x i x i+2 ) • 

i=l ^ i=0 ' 

(5) 

Using standard techniques for diagonalizing such 
Hamiltonians [25] (and the trick of adding an extra qubit 
to make the system have quadratic fermion operators, 
see Appendix A), this Hamiltonian can be shown to be 
equivalent to a system of non-interacting fermions with 

Ms) 2 = 4 (l - 2s(l - s) [l - cos(^)] ) (6) 

where r]\, is the creation operator for the hth. fermion 
(see Appendix A) and I = [|J or / = [^^J depend- 
ing on whether one is considering the even or odd chain. 
Note that there is no k = energy level. From this 
equation one sees that the minimum energy gap occurs 
at s = \ where 0^(1/2) = 2 cos( 2 (f+i) ) and is of order 

O(-)- Thus, since the minimum gap scales inversely as a 
polynomial in the length of this one-dimensional system, 
we see that if we turn the field on in a time polynomial 
in this length, then with high probability the quantum 
information will propagate from one end of the system to 
the other end — with a possible Hadamard gate applied, 
depending on whether n is odd or even. 
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B. Single Qubit Quantum Circuits 

Having demonstrated rigorously that the scaling for a 
quantum wire has a gap which scales inversely with the 
length of the wire, we now examine what happens when 
we apply other single qubit gates by using twisted Hamil- 
tonians with varying 0's. Here we give strong numerical 
evidence that the gap scales inversely as a polynomial 
through the use of a matrix product state algorithm. To 
simplify our study of the one dimensional twisted Hamil- 
tonian, it is convenient to note that this model has a 
duality. In particular, for a twisted Hamiltonian with 
angles Oi where 1 < i < n — 2 and Oi is the angle associ- 
ated with the (i + l)st qubit, it can be shown that if the 
angles obey the symmetry condition Oi = 6 n -i-i for all i 
and additionally there is only one minimum in this model 
(which is true numerically for small systems), then the 
minimum energy gap occurs at the midpoint, s = 1/2. 

To see this, we proceed as follows. First we define a 
shorthand for the operators which are rotated combina- 
tions of X and Y Pauli operators, namely, 

[e k } i = cos(e k )[X} i + S m(6 k )[Y} i . (7) 

Now consider the twisted one-dimensional cluster Hamil- 
tonian with a transverse field, 

H(a) = - (1 - s)(j^]iNi+i[^+2 + [Z] n -i[X))j 
-s J2lX]i- (8) 

i=l 

Define the following unitary operator: 

U = f[[X]^j (j[[C z ] iA+1 j QjWflOli+iJ 

(9) 

where R(0) = exp(— iZO/2), [Cz]ij is the controlled- 
phase gate acting between qubits i and j, and S is the 
gate which inverts the chain about its middle, swaping 
the 1st and nth, 2nd and (n— l)st, etc. qubits. Then one 
can check that 

xn-2 

UH(s)rf =~s(J2 i Z U n-i-l}i+l [Z] i+ 2 

^2=1 

)n — 1 
-(1 (10) 

Thus, in the case that Oi = # n _^_i for 1 < i < n — 2, H(s) 
has the same spectrum as H(l — s) and in particular if 
there is only one quantum phase transition (meaning, 
only one value of s for which the spectral gap above the 
degenerate ground space collapses to zero in the thermo- 
dynamic limit), then it occurs when 8 = 1 — s or s = \. 
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FIG. 5. Energy gap versus number of qubits for a one- 
dimensional twisted Hamiltonian with random angles. The 
duality condition has been enforced so that the minimum is 
at s = | (see text.) 200 samples were taken with uniformly 
distributed random angles. We also plot the smallest and 
largest energy gaps that were obtained for a given number of 
qubits (dotted lines) as well as the case with no rotated an- 
gles (which we showed could be computed exactly) (non-bold 
line.) The data is well fit by (4.97 ± 0.02)n _1 . 



Given this, we used the technique of matrix product 
states [26-28] to calculate the energy gap for systems 
having up to 200 qubits with angles 0i chosen uniformly 
from [0, 2ir) satisfying the aforementioned symmetry con- 
dition; for each system size that was examined, we per- 
formed 200 simulation runs, each with a different choice 
of angles. Fig. 5 presents the resulting data, from which 
we find that the energy gap data is well fit by a function 
which scales as 0(n _1 ), thus indicating that the twisting 
of the Hamiltonian does not quantitatively change the 
scaling of the energy gap with the length of the single 
qubit circuit being implemented. 



C. Mult-qubit Quantum Circuits 

Finally we turn to the question of the energy gap in 
the case where the circuits involve two or more qubits. 
Here we note that if we take a square lattice with ap- 
propriate boundary conditions, then with the untwisted 
Hamiltonian model is equivalent to a quantum compass 
model restricted to a certain symmetry sector. Under 
this symmetry restriction, the minimal gap occurs at 
s = |, where the compass model has equally competing 
interactions. Numerical evidence from exact diagonaliza- 
tion on a square lattice for this model presented in [29] 
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demonstrates that the minimal energy gap for this model 
with the symmetry restriction scales as 1/n where n is 
the size of (i.e., number of qubits in) the square lattice. 
Thus the real question is what happens to the energy gap 
for more general circuits. 

We consider a model on a square lattice with the full 
twisted cluster state Hamiltonian turned on interpolating 
to a Hamiltonian with applied fields on across the entire 
device. Label the qubits on the square grid by (z, j). The 
Hamiltonian we consider is 



H =(1 - s)H { 



c 



(11) 



where L is the length of the lattice and we define [P]^j\s 
with values of i or j lying outside the lattice as identity 
operators to deal with boundary terms. Here, He is the 
twisted cluster Hamiltonian of Eq. (1) for a square lattice 
with some particular choice of angles 6uj\ for each vertex 
(i, j) in the lattice. We suppress the dependence on the 
0(ij) to avoid notational clutter. 

If we apply a unitary operator U x consisting of a 
controlled-phase gate between all ^-coordinate neighbors, 
[Cz](i,j),(i+i,j), to the Hamiltonian in Eq. (11), then the 
cluster state Hamiltonian is turned into a y-direction 
striped twisted Hamiltonian and the applied field turns 
into a x-direction striped untwisted Hamiltonian. Explic- 
itly, this transformed Hamiltonian (which will have the 
same spectrum) is given by 

L 

+ s i Z ](ij-l)[ X ](i,j)[ Z ](iJ+l) ■ (12) 

If one instead applies U y , consisting of controlled-phase 
gates between all ^-coordinate neighbors, [Cz](i,j),(i,j+i) 
to the Hamiltonian in Eq. (11), then we see that this 
simply swaps the direction of the stripes: 
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FIG. 6. Energy gap versus number of qubits for the two- 
dimensional twisted Hamiltonian with random angles and the 
duality condition enforced (see text.) The initial Hamilto- 
nian was a rotated cluster Hamiltonian for a square lattice 
(Eq. (1)) and the final Hamiltonian was an X field over the 
entire device. All data were obtained using exact diagonal- 
ization with accuracy 10 -8 , with the exception of the random 
angle data for n — 25 which were obtained using a one di- 
mensional matrix product algorithm wrapped onto the two- 
dimensional lattice and has accuracy 10 -3 . The data is well fit 
by (1.63 =b 0.12)n _1 . We also show the untwisted energy gap 
for comparison. The untwisted energy gap does not appear 
to be a power law because there is an even-odd lattice effect: 
examination of non-square lattices gives us more confidence 
that the untwisted case is indeed a power law. 



III. FAULT-TOLERANCE OF THE ADIABATIC 
QUANTUM TRANSISTOR 



U y HU l = ~ Yl ( x ~ s )\. (ij)hj)l Z hj-i)l Z hd+V 

+ s [ Z ](i-lj)[ X ](ij)[ Z ](i+l,j) ■ ( 13 ) 

If one swaps the qubits about the line x = y in 
Eq. (13), we see that we obtain the striped Hamiltonian 
in Eq. (12), but with 0^j\ replaced by and s re- 

placed by 1 — s. Thus if we enforce O^j^ = we 
obtain a duality, and as in the previous subsection there 
will be a minimum at s = \. 

For systems which obey the duality condition 0(ij) = 
we have investigated the size of the energy gap using 
exact diagonalization as well as using a matrix product 
state approach. The results of these simulations are plot- 
ted in Fig 6. While we can only obtain weaker evidence 
of an inverse polynomial for this two-dimensional system, 
the data are at least consistent with this hypothesis. 



Having given strong evidence that our model can effi- 
ciently produce a desired quantum circuit in a case where 
all evolutions are error free, we now turn to the question 
of how this model will behave in realistic settings where 
the system is coupled to an environment. Here we argue 
that if one uses standard circuit constructions for fault- 
tolerant QC [3-5], that this will be good enough to ensure 
robust quantum computation in the model. In particular 
we can consider a model in which each modular adiabatic 
quantum transistor is configured to perform an encoded 
quantum gate, state preparation, or measurement from 
a fault-tolerant quantum circuit construction, including 
the error correcting step. We will argue that if one does 
this, then the standard analysis of the success probabil- 
ity of such constructions carries over to our model. Note 
that fault-tolerant QC requires expunging entropy (usu- 
ally via measurement), but that this can always be placed 
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at the end of the fault-tolerant block [30]. Thus we en- 
vision here a model in which the final measurements in 
an error correcting circuits are implemented after each 
evolution of a quantum transistor. This can be done by 
either directly measuring the relevant qubits, or via an 
adiabatic amplifier which dissipates energy by the natu- 
ral relaxation of the system, as described in Sec. VA. 

First note the following good news about using each 
quantum transistor as a fault-tolerant circuit gadget. If 
one is executing a quantum algorithm with L gates, then 
each of these gates needs to be executed with accuracy 
e = eo£ -1 . To do this using fault-tolerant gates below 
the fault-tolerant threshold requires [3-5] that the gad- 
gets have circuits of size 0(log c e _1 ) for a constant c. 
Thus if the energy gap in an adiabatic gate shrinks as 
an inverse polynomial, the energy gap in the individual 
fault-tolerant gadget will shrink as 0((log fc L) _1 ) for a 
constant k. Because each individual transistor executes 
a single encoded gate, this means that each adiabatic 
evolution need only be polylogarithmically longer when 
executing larger and larger quantum algorithms. 



A. An Unrealistic But Illuminating Error Model 

Now turn to the question of whether fault-tolerant gad- 
gets will act to overcome errors in the adiabatic transis- 
tor. Begin first by considering the following simplified 
error model: assume that at the beginning of the com- 
putation the system is not in the ground state of He, 
but rather is in one of the excited states of He- Note 
that since the terms of He commute (as expressed in 
Eq. 1), excited states of He can be labeled by the list of 
eigenvalues of each of the terms, where the eigenvalue of 
each term is either +1 or —1. If one carries through 
the analysis described in [16] regarding the computa- 
tion that is performed if one starts in such an excited 
state, then one sees that for each of the terms in He 
whose initial eigenvalue is +1 there is a corresponding 
Pauli error in the circuit — in other words, errors that 
take the form of starting in an excited state of He map 
directly onto errors in the circuit model. For example, 
consider the one-dimensional case with no rotated an- 
gles; if the initial state is in the +1 eigenspace of some 
term — [Z]i_i[X] i[Z]i+i in the middle of the chain, then 
this corresponds to a Pauli Z error in the circuit model 
of this circuit. 

Furthermore, note that in this model errors which are 
local on the physical qubits of He map to errors that are 
local in the circuit since all terms in He are localized on a 
few qubits — so for example, applying [X]i to the system 
flips only two eigenvalues of He in the untwisted case. By 
the linearity of quantum mechanics, we can extend this 
argument to initial preparations which are mixtures and 
coherent superpositions of such errors. Thus if we take 
our initial system and expose it to an environment for an 
amount of time proportional to the size of the quantum 
circuit we are about to enact, then subsequent perfect 



adiabatic evolution will produce errors in the quantum 
circuit model which are commensurate with a local inde- 
pendent error model. (Note that this model does have 
temporally correlated errors over a few qubits, but this 
does not affect the existence of fault-tolerant methods for 
overcoming these errors [3].) 

B. Errors On A Quantum Wire 

Of course, this model in which errors appear only at the 
beginning of the computation is a fiction. More generally 
we may consider a model in which errors occur through- 
out the adiabatic evolution. Consider such an evolution 
for the single quantum wire with an untwisted Hamilto- 
nian. As we detailed above, this model can be mapped 
into two transverse-field Ising models. Any linear opera- 
tor, and hence any error, on this model can be expressed 
as one of three types of errors: errors that change the 
energy of the system, errors that act on the degeneracy 
(corresponding to the encoded qubit), and a combination 
of these two errors. We will show that each of the first two 
types of errors can be mapped to an independent error 
model on the quantum circuit corresponding to the quan- 
tum wire in an independent manner, thus taking care of 
the third type of error. Finally note that we use the term 
error in a quantum error correcting sense to denote an 
error operator, even when this arises from coherent error 
sources. For example our results deal with the errors in 
the Hamiltonian description resulting from a perturbing 
interaction on the system H' = H + XV for small A as 
well as for couplings of our system to an environment, 
though the argument for these errors is slightly different 
(see Sec. IIIB3 for a discussion of this point.) 

1. Errors that change the energy of the system 

We first consider errors that change the energy of the 
system. The energy spectrum calculated in Sec. II A has 
low-lying excited states and thus we can reasonably ex- 
pect excitations into these levels. Assuming detailed bal- 
ance in our error rates (which tends to be the case for 
weak coupling errors that change the energy of a sys- 
tem [31]) we consider a model in which the rate of errors 

that change energy by AE scales as p^e k B T where po is 
a bare error rate, ks is the Boltzmann constant, and T 
is the temperature of the environment. 

Consider working at a temperature much less than the 
bare energy gap of the system at s = 0, which is A = 2. 
In this case the only fermionic excitations that will occur 
with high probability are those for which u(s) < ksT. 
This condition implies that only a constant fraction of 
fermions will be created in such an error model through- 
out the adiabatic evolution. Furthermore, note that be- 
cause the adiabatic evolution preserves the energy levels, 
it will therefore proceed to drag the system along the 
excited energy level. We see therefore that errors that 
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FIG. 7. When one exactly solves the one- dimensional wire, 
the result is that the wire behaves like a collection of free 
fermions. Here we show the energies for creating each of these 
fermions in a chain of length 80 (thus all the energy levels are 
sums of these energies) as a function of the scaled time s. 
See Eq. (5) in Appendix A. For high energy excitations, if 
the environment is cool enough, then these energy levels are 
essentially never excited (solid lines.) At energies comparable 
to the temperature of the environment, errors will no longer 
be suppressed by a Boltzmann factor (dashed lines.) At the 
lowest energies, the interaction strength between the environ- 
ment and the system is small enough that the error model 
will not be governed by detailed balance. However, assuming 
that these energy levels are changed by any process whatso- 
ever is no worse than our assumption that all errors below 
the temperature cutoff lead to real errors. 



excite a constant fraction of fermions can be mapped 
back to excitations at the beginning of the adiabatic evo- 
lution, which we have already shown are not a problem 
for standard quantum error correction procedures. 

Finally we note a technical condition: we have assumed 
a detailed balance error model, but for the lowest lying 
energy levels, where the strength of the interaction be- 
tween the system and environment is greater than the en- 
ergy gap, this condition might not hold. However, these 
energy levels are, under our temperature assumption, al- 
ready assumed to be in error, and therefore covered by 
our argument (i.e., they are at most a constant fraction of 
the error, the fraction being directly related to the ratio 
of the energy gap and the perturbation energy strength.) 
See Fig. 7 for a graphic depicting the relevant energy 
levels and the partition of these errors into different cat- 
egories in our argument. 

To summarize: under the reasonable assumption that 
high energy excitations are suppressed by a Boltzmann 
factor, we see that errors in the quantum wire that change 
the energy of the system result in independent errors on 
the quantum circuit being enacted. 



2. Errors on encoded quantum information 

The situation for errors that do not change energy but 
instead act on the degeneracy of the quantum wire is 
slightly more complex. To see what happens, we con- 
sider the two logical operators encoded in the degenerate 
ground state of the quantum wire, which for a chain of 
odd length n are given by 

n-\-l re.— 1 

2 2 

X n =l[[X] 2i _ 1 and Z n = [Z] n J[[X] 2i . (14) 

The first n — 1 operators are the encoded operators of the 
code described in Sec. II A. Note that these two logical 
operators commute with the entire Hamiltonian during 
the adiabatic evolution. At the beginning (end) of the 
computation this information is localized onto the start 
(end) of the quantum wire. Thus during both of these 
times the quantum information in the wire is susceptible 
to decoherence from local errors. For example undesired 
or imprecisely controlled terms in our initial Hamiltonian 
will be able to act non-trivially on the initial quantum 
information. However, during the middle of the compu- 
tation, we will argue that the information in the degen- 
eracy is protected, and furthermore that we can adjust 
the adiabatic schedule in such a way that the system only 
spends a constant amount of time during which the quan- 
tum information is exposed to local decoherence at the 
beginning and end of the adiabatic evolution. 

In particular consider the logical X n . By expanding 
this logical operator in terms of the underlying physical 
operators, we see that this error requires a combination 
of 0(n) local [X]i errors on odd qubits. When we map 
these errors over to the transverse Ising model, these are 
errors on the Ising model, but are now two-qubit terms 
like XiXi+i for i > 2 or X\. Individually these errors can 
change the energy of the system (and we have previously 
argued that such errors can be dealt with), but they can 
also have an effect which keeps the energy constant. To 
evaluate this effect, for example, for information at the 
end of the computation, we can evaluate the portion of 
the error amplitude which preserves the vacuum of the 
transverse Ising model: e = ^^(^1X^X^+1 \v) where \v) 
is the vacuum (the single X\ error can be thought of as 
a two qubit X error on the fictitious extra qubit.) 

For the transverse Ising model with periodic bound- 
ary conditions in the thermodynamic limit, the nearest 
neighbor XX correlations are given by [32] 

(v\[X] i [X] i+1 \v) = - [ IT^cos^+ajdfc, (15) 
71 Jo 

where 

ft(Jfe) = v / l + a 2 + 2acosfc (16) 
and, to relate back to the quantum wire model, a = 
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This can be transformed into 



of 



(v\[XUX] i+1 \v) 



1 

7TS 



(2s-l)K[4s(l-s)] 
1 

ITS 



-E[4s(l-s)] 



(17) 



where K and £" are the complete elliptic integrals of the 
first and second kind respectively. We are interested 
in the amplitude for the product of n of these correla- 
tion functions, (v\ [X]i[X]i+i \v) n ^ 2 . Series expanding this 
about s = 1, we find 

(v\[XUX} i+1 |t;)"/ 2 = 1 - In(* - l) 2 + \n{s - if 



2n 2 -53n\. 4 
+ 0((s-lf). 



(18) 



If < 0(n"2- e ) for e > 0, we see that this function 

vanishes as n -+ oo. 

Thus using correlation functions for the relevant de- 
generacy preserving operation we see that the amplitude 
for the degeneracy-preserving error is exponentially small 
(as a function of n) except for s near s = 1 in a window of 
size 0(n~ 2 _e ) for e > 0. In and of itself this would imply 
that the model is in trouble: with a linear interpolation 
for this would imply spending 0(y/n) time in this region 
during which the quantum information can be decohered. 
However, note that the energy gap for this model is large 
during the beginning and end of the adiabatic path. Thus 
instead of using a linear interpolation as in Eq. (2) one 
can use an interpolation which spends only a constant 
amount of time in the window where these errors can 
affect the logical quantum information. 

In particular if we maintain the adiabicity condition for 
each infinitesimal evolution of the system we can obtain a 
schedule for interpolating our Hamiltonian which spends 
significantly more time where the energy gap is smaller. 
In particular, consider the adiabatic condition for the 
one-dimensional quantum wire. As we have shown above, 
the energy spectrum for a wire of length 21 is given by 
Eq. (A22). To maintain the adiabatic condition for the 
lowest energy level infinitesimally we must satisfy [33] 



ds 
dt 



dt 



ds 



(19) 



where e is the accuracy that we desire for our adiabatic 
evolution relative to the ideal evolution. Performing this 
integration and enforcing t = at s = yields the solu- 
tion 



1 

t =— CSC 
4e 



(— 

irl 



2(1 + 1) 



tan 



(2s — 1) cot (2(1+1)) 



(20) 



Inverting this equation for s yields a schedule for adjust- 
ing s as a function of time which runs for a total time 



T = 



ttI 



4e(Z + 1) 



CSC 



(+) 



(21) 



which scales linearly as I for large I (thus there is no 
Grover-type speedup using this schedule.) The amount of 
time during which this schedule spends with 1 — 77 < s < 1 



2e 



2(Z + 1) 
ttI 



tan 



(2r] - 1) cot 1 



2(2+1) 



(22) 



For the case relevant to the degenerate error model where 
r] = l~ x l 2 this is bounded by K Thus this schedule 
spends only a constant amount of time in the region 
where the logical output quantum information is suscep- 
tible to noise. 

In conclusion, for the single qubit quantum wire case, 
we see that the quantum information in the degeneracy 
is vulnerable to degeneracy-preserving errors only at the 
beginning or end of the computation, and in a setting 
which is no worse than the standard circuit model with 
independent errors. We can adjust the adiabatic schedule 
to not spend much time in these regions. Really, this is no 
different than how we would want to operate a quantum 
computer with or without adiabatic quantum transistors: 
to keep the effective error-corrected error rate down, one 
needs to spend a minimal amount of time letting the 
quantum information simply decohere. In our model this 
means that we need to move faster than linearly out of 
the region where the encoded quantum information is ex- 
posed to errors. We note that the above analysis works 
for both open system errors and also for (small) devia- 
tions in the Hamiltonian we implement, i.e. perturbative 
system errors. 



3. Static Errors 

In the previous subsection we have argued that an adi- 
abatic quantum transistor configured to act as a quan- 
tum wire would behave, given an appropriate adiabatic 
schedule, as a quantum circuit with an independent error 
model. In this argument we have explicitly used an ar- 
gument that relied on detailed balance. Thus one might 
wonder how the argument works when there is no envi- 
ronment and the errors occur solely from a static per- 
turbation on the wire. Here we argue that these errors 
also give rise to an independent error model, at least for 
a quantum wire with no gates. 

In particular consider a model in which a static term 
is added to our adiabatic evolution H'(s) = H(s) + XV 
for the case of a quantum wire. We will assume that V is 
a sum of local operators (for example, one or two qubit 
operators) and that the strength of the interaction is per- 
turbative (A <C 1, since we have fixed to unity the energy 
scale used in the twisted cluster-state Hamiltonian and 
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the final applied field.) Now first note that the contri- 
bution of V to splitting the degenerate ground state of 
the wire will be exponentially suppressed except at the 
beginning and end of the adiabatic evolution, as we have 
argued in the main text. The reason for this is exactly 
the same: local operators do not act on this space except 
when multiplied together (except at the beginning and 
end of the adiabatic evolution.) Thus a term like XV will 
act like a single-qubit error term only for the time we 
spend with the information localized on either end of the 
wire. Using the scheduling trick we have described pre- 
viously, this means that it acts like an independent error 
(from a static term) at the beginning and end of the com- 
putation represented in the wire. Thus we see that the 
static perturbing term will not act on the degeneracy in 
a manner that is inconsistent with an independent error 
model. 

Moreover, many types of errors will hardly effect the 
degeneracy at all since it is known [34] that the cluster 
state has symmetry-protected topological order [35, 36] 
and the degeneracy and spectral gap for such systems is 
known to be stable [37, 38]. (See also Refs. [39, 40] for 
other schemes for quantum computation in a symmetry- 
protected topologically ordered phase.) 

Next we turn to the case where there is a static per- 
turbation, but where we now concern ourselves with the 
effect this has coming from portions of V that do not act 
on the degeneracy. These have the effect of changing the 
energy levels, both in terms of the eigenvalues and the 
eigenvectors of the system. 

Consider first the effect of the perturbation on the 
eigenvectors of the system. We begin by recalling the ar- 
gument that shows that the quantum wire Hamiltonian 
correctly transmits information in the absence of a static 
perturbation. Recall that the quantum wire Hamiltonian 
is given by 

H(t)=-f(s)(Aj2[Z] i [X} i+1 [Z] i+2 -A[Z} n _ 1 [X] n ^j 

n-1 

- g(s)Aj2[X] i . (23) 

Initially we can define the logical operators X = [X]i [Z] 2 
and 2 — \Z\\. Initially the system is in the +1 eigenstate 
of the vertex operators [Z]i[X]i+i[Z]i+2 (1 < i* < n — 2) 
and [Z] n _i[X] n . By suitable multiplication of the en- 
coded logical operators by these vertex operators, we 
can express the logical operators as a pattern of [X]i 
operators on qubits, for i < n, and a Pauli operator 

on the last qubit: X = (U7=i,i odd!*]*) l X U, Z = 

(jXi=ii evenM^) [^]n? where we have assumed for sim- 
plicity that n is odd. At the end of the evolution we will 
be in the +1 eigenstate of all of the [X]i, 1 < i < n — 1. 
This then implies that the information will be correctly 
transmitted down the quantum wire since over this con- 
dition the logical operators are just bare Pauli operators 



on the last qubit. 

Suppose, on the other hand, that we were initially in a 
state in which we were in the —1 eigenvalue eigenspace of 
a single vertex operator. Then the above argument about 
the logical operators could still be applied, but the effect 
would now be that one of the logical operators might ac- 
quire a phase. For example it might end up that X = 

-(nKioddMi)[*]», 2 = (TCm events) 

This then represents the quantum wire being faithful up 
to a Z error. Generalizing, we see that initial states that 
are not in the all +1 eigenstates translate into single qubit 
errors on the wire (so that multiple vertex operators be- 
ing in —1 eigenstates represent multiple errors.) A similar 
argument applies to the end of the evolution: if we are 
are in —1 eigenstates of some final [X]i operator, this will 
correspond to a single-qubit error on the wire. Note that 
it is possible for the errors in our model to cancel out. 

How does the static perturbation changing the eigen- 
values affect this argument? Consider first the effect of a 
perturbation on the initial ground state and expanded its 
effect using perturbation theory. In particular if \ip g (0)) 
is the initial ground state of the unperturbed system 
|V^(0)) is the initial ground state of the perturbed sys- 
tem, then the first order correction will be 

W(0))H^)) + aEM»W0)) 

fc ^ E g -E k 
+ 0(A 2 ) (24) 

where 1^(0)) are the excited energy states of the unper- 
turbed system. Now notice that if V is a sum of local 
terms, each local term will act to change the eigenvalues 
of at most three vertex operators. Thus the sum above 
will not be over all excited states, but will only included 
excited states with at most three vertex operator eigen- 
values flipped. Each of these terms represents errors that 
are localized in space-time on the quantum circuit version 
of the quantum wire, as per our argument in the previous 
subsection. If we were to end up in the +1 eigenstate of 
the final Hamiltonian's [X]i terms, these would then each 
represent a wire in which at most three single qubit er- 
rors occurred during the evolution. That is, the portion 
of the wave-function arising to first order in perturbation 
theory, if dragged to the +1 eigenstates of the [X]i oper- 
ators, acts as if one had (nearly) independently erred the 
quantum information transmitted through the quantum 
circuit. A similar argument can be made for the effect 
of being in the wrong final state. Thus we see that the 
effect of being in the wrong initial and final eigenspace 
is, for a weak local perturbation, equivalent to an error 
that is local in space and time in the quantum circuit. 

Having shown that the effects of changing the eigenvec- 
tors of the ground state can be modeled as an indepen- 
dent error model on the quantum wire, we now turn to 
the effect of static perturbations on the eigenvalues. The 
worry here is as follows: the effect of the perturbation 
may cause the energy gap to close and therefore the adi- 
abatic evolution will actually cause the evolution to not 
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FIG. 8. Static perturbations. Consider a quantum wire with 
repeated Hadamard gates, but where the individual interac- 
tions in, for example Eq. (23) are randomly perturbed by a 
random amount in the range of— A to A. This plot is done for a 
wire of length 20. This model can be mapped to the transverse 
Ising model and above we plot the energy of the free fermions 
that occur in this model. One sees that while the energy cost 
for each fermion with the perturbed field is modified qualita- 
tively the cost for each individual fermion remains about the 
same. This implies that the model will behave similarly to the 
unperturbed model: in the absence of other errors the worst 
this can do is cause a level crossing near s = 1/2 which will 
at most introduce a constant number of errors in the model. 
Further the effects of other errors occurring on this perturbed 
model will behave similarly to how we argue they will act in 
the main text. In particular there will be a constant fraction 
of errors occurring when the system is coupled to a thermal 
bath and obeys detailed balance. 



preserve the ground state (subspace). This will certainly 
happen in the system, but the real question is whether 
this looks like an independent error model or not. 

To see that this is unlikely to be a problem, consider 
a model in which each of the individual terms in the 
Hamiltonian of the quantum wire, Eq. (23) are randomly 
perturbed by ±A. When we convert this model over into 
the free fermion model as discussed in previous sections, 
we see that this will result in A and B matrices from 
Eq. (A3) that are like the band diagonal matrices in 
Eqs. (A5) and (A6), but now with random diagonal and 
off-diagonal elements added to these matrices. One can 
then numerically diagonalize the equation for the free 



fermion energies. In Fig. 8, we show both the unper- 
turbed free fermion energies and also the perturbed en- 
ergies for a chain of length 20 and a perturbation strength 
of A = 0.1. From this one see that while the perturbation 
causes the energies associated with the free fermions to 
change by ±A, the spectrum is qualitatively the same. 
In particular we can now imagine sweeping adiabatically 
from s = to s = 1. The energy levels that will cause 
problems are then the free fermions with excitation en- 
ergies below A. At all stages where the fermion energies 
are greater than A there is a cost to create such a fermion 
which is greater than A, and the perturbing potential V 
will not be strong enough to achieve this transition. Since 
there are only a constant fraction of the energy levels that 
cross below the A creation energy line (similar to the main 
text Fig. 7) this means that the possible errors created 
are a constant fraction of those that could occur on the 
quantum wire. This, then, is nothing more than an in- 
dependent error model on the system. Thus, at least for 
the kind of error model we have assumed, static pertur- 
bations result in an independent error model. 

We have argued (but not rigorously proven) in this 
subsection that the effect of a static perturbation on a 
quantum wire is no different from an independent error 
model on the quantum circuit corresponding to this wire. 



4 . Summary 

In this section we have argued that errors on quantum 
adiabatic transistors configured to act as a quantum wire 
act similar to errors within the standard quantum circuit 
model with independent (or slightly time correlated) er- 
rors. Errors that change energy were argued to only act 
as a constant fraction of independent errors in the quan- 
tum circuit model. Errors that preserved the degener- 
acy could be severe, but only during the initial and final 
durations of the adiabatic evolution when the quantum 
information is localized on either end of the quantum 
wire. By spending less time in the adiabatic evolution 
for this time we showed that these errors act as constant 
rate errors at the beginning and end of the computation. 
Finally we analyzed how static (not varying in time) per- 
turbative interactions would modify these arguments and 
gave evidence that these errors behave in a similar model. 
Thus adiabatic quantum transistors configured to act as 
a quantum wire seems to have errors that act as indepen- 
dent single qubit errors in the standard quantum circuit 
model. The case for a full quantum circuit with multi- 
ple qubits is more difficult to analytically or numerically 
analyze. Some of our arguments, for instance the effect 
of local perturbations on the initial twisted cluster state 
Hamiltonian being local in the space-time of the quantum 
circuit enacted by the adiabatic evolution, carry over di- 
rectly to this more complicated case. The the full general 
case is an open question for future work. 
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IV. PERTURBATION THEORY GADGETS 

Finally let us mention issues arising from the use of 
perturbation theory gadgets in implementing the many- 
qubit interactions in He- These gadgets produce an ef- 
fective, low-energy many-qubit interaction via strong and 
weak single- and two-qubit interactions, the strength of 
these interactions being related to the ratio of the weak 
and strong interactions (hence the name.) If one wishes 
to only apply fields (single-qubit interactions) to the sys- 
tem, it appears that one cannot use Bartlett-Rudolph 
gadgets [20], but instead must use gadgets that do not 
use encoded quantum information [21] (this is because 
the former gadgets have encoded quantum information 
that is not acted upon by single qubit operations.) At 
temperatures much lower than the energy gap in the an- 
cillas of these constructions, one can effectively treat the 
errors arising from these gadgets as local errors in our 
above model. In particular, as in prior work [15, 20, 21], 
we can treat the imperfections arising from the use of 
perturbation gadgets as terms that cause errors on our 
system. Thus the only major quantitative change result- 
ing from using these gadgets is that the strengths of the 
engineered interactions will be lower due to the pertur- 
bative nature of these gadgets. The best possible choice 
of gadget will depend on the physical system being used 
and on the relevant noise sources that are present. 

The main question which arises from the use of per- 
turbation theory gadgets is whether they will cause extra 
problems when the quantum transistor is configured to 
carry out a fault-tolerant quantum circuit element. In 
the supplementary section of Ref. [15], we analyzed the 
effects of perturbation theory gadgets for implementing a 
simple adiabatic scheme that performs a two-qubit gate. 
The main result of this analysis was that the use of per- 
turbation theory gadgets to achieve these many-body in- 
teractions does not substantially damage the usefulness 
of these gadgets in the adiabatic constructions. However, 
the use of gadgets has two drawbacks: (1) the eigenvec- 
tors are slightly changed due to the inexact nature of 
the effective many-qubit interaction, and (2) the eigen- 
values are also slightly changed and thus the gap is made 
slightly smaller. For the case analyzed in the paper, there 
were two energy scales, A and A. We can choose units 
where A = 1. The effective three-qubit interaction in the 
scheme was of strength 0(A 2 ), and thus the gap in this 
construction is of this order. The imperfections in the 
gadget construction led to corrections to the energy of or- 
der 0(A 3 ). Furthermore, the eigenvectors of the ground 
state were corrected in amplitude to order 0(A 2 ). The 
main conclusion of this calculation was that the effects 
of the perturbation theory gadgets was small in the per- 
turbation parameter A. 

What can we say about the similar calculation for the 
adiabatic quantum transistor? The first question is how 
the gadgets will affect the initial and final eigenvectors 
of the system. Bartlett and Rudolph [20] considered ex- 
actly this problem for their encoded gadget scheme, and 



showed that the new eigenvectors behave as if there was 
an independent error model acting on the state when 
it was used for measurement-based quantum computing. 
One can perform exactly the same sort of calculation 
for the perturbation theory gadgets that use ancillas as 
mediators and the conclusion is exactly the same: the 
lowest-order errors created using these gadgets change 
the state to act as if it has been independently erred. 
The result is not dissimilar to that of our Eq. (24), but 
now for the more complicated quantum transistor with 
a quantum circuit. The other question is what effect 
the gadgets have on the energy gap. We note here that 
the corrections to the effective interaction for the pertur- 
bation theory gadgets can be thought of as extra static 
Hamiltonian terms that are weaker by the perturbative 
factor used to create the gadgets, hence these errors are 
covered by the argument in the previous section. 

The use of perturbation theory gadgets in building adi- 
abatic quantum transistors is clearly one of the largest 
drawbacks in our scheme and one of the reasons why 
we consider this work as an outline for how a adiabatic 
quantum transistors will work. Recent work using, for 
example, more physically realistic AKLT states [39-43] 
however indicates that these gadgets may not be a nec- 
essary component of a quantum transistor construction. 



V. BUILDING BLOCKS FOR ADIABATIC 
QUANTUM TRANSISTOR ARCHITECTURES 

Finally let us detail some building blocks for adiabatic 
quantum transistors that will be useful for building a 
larger quantum computer architecture that closely mim- 
ics today's modern clocked classical computer architec- 
tures. These include constructions for a system that per- 
forms measurement, and a system that can be used to 
take classical information and use it to route quantum 
information in a device built with an adiabatic quantum 
transistor. The former is essential for fault-tolerant con- 
structions, which must purge the entropy of quantum er- 
rors, and the later is important for building architectures 
that are modular and synchronous. 



A. Adiabatic measurement amplifier 

Here we describe the details of an adiabatic measure- 
ment amplifier. This is a way to take quantum infor- 
mation encoded on a single quantum wire (in the circuit 
being adiabatically simulated) and copy it to multiple 
qubits in a fixed basis and then read out this informa- 
tion. Consider a tree with alternating internal levels of 
degree 2 and 3 as shown in Fig. 9. Call the degree one ex- 
ternal node that is connected to an internal degree three 
node the root of the tree, r, and its one child, r' and let 
the other degree 1 vertices be called leaves. Then ini- 
tially information will have logical operators Z = Z r and 
X = X r Z r > . The operation of the amplifier is as follows. 
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FIG. 9. An adiabatic measurement amplifier. A binary tree 
with alternating degree 2 and and 3 internal vertices (as 
shown) can be used to take information encoded into the X 
eigenstates at the root into a code across the leaves by turn- 
ing on an X field across the shaded qubits. This code has a 
stabilizer given by ZiZj terms on leaves, i.e. it is the repeti- 
tion code. The logical information initially at the root in the 
X eigenstate is in this process copied into the Z eigenstates 
of the leaves. Thus If one turns on a Zi interaction, the 
information initially in the X eigenstate will be separated in 
energy by twice the number of leaves. The natural relaxation 
of this system to it's ground state will result in 0(# leaves) re- 
laxation events, which can be individually (noisily) observed, 
hence amplifying a measurement in the X eigenvector basis. 



First, single-qubit X terms are turned on while turning 
off all of the graph stabilizers. Then a Zi term is 
turned on across all of the leaves. One then waits for 
observed decay events from the leaves: no decays indi- 
cate the system was in the —1 eigenstate of X initially, 
0(#leaves) decays indicate the system was in the +1 
eigenstate of X initially. The utility of this devices is 
that it turns information encoded into one qubit into in- 
formation copied (in a particular basis) to many qubits, 
thus facilitating the measurement of this information. 

To see how this works, we note that while we defined 
the information as initially encoded into Z = Z r and 
X — X r Z r t , we could have chosen equivalent operators 
by multiplying by graph stabilizers. Thus if we multiply 
the Z operator by graph stabilizers from vertices of de- 
gree 3 and the graph stabilizer from the leaves, we obtain 
an operator that has X or I on the non-leaf vertices, and 
is X on all of the leaves. Similarly, if we multiply a X 
operator by a graph stabilizer from vertices of degree 2 
on any simple path from a root to a leaf we will obtain 
X operators on the non-leaves, and a Z on the leaf. This 
later fact is independent of which path to the leaf one 
choses. This is a symptom of the fact that the informa- 
tion is encoded into an error correcting code. Indeed it 
is not hard to see that by multiplying graph stabilizer 
elements one can obtain ZiZj terms that act only on i 
and j being leaves of the graph. These are the stabilizer 
elements for the classical repetition code: \b) —> \b)® k . 
When we turn on the X terms on the non-leaves we will 
thus have taken the information at the root (and its child) 



and encoded it into this code (note that the role of X and 
Z are reversed in this process due to an extra Hadamard 
at the beginning of this process). Stabilizer code argu- 
ments [15, 16] make this statement rigorous: we promote 
the graph stabilizer to logical Zs and then the local X 
terms become products of logical Xs. The information 
as described by the encoding above is untouched by this 
process and thus the adiabatic evolution described above 
does not affect this information. 

An equivalent way to derive this graph is to take 
the circuit for copying information recursively using 
controlled- not s and prepared ancilla states, convert it to 
a graph state, and then simplify this state by noting, for 
instance, that two nodes on a line can be eliminated from 
a graph because they correspond to H 2 = J, where H is 
the Hadamard gate. 

B. Adiabatic Router 

It is convenient to design a device that can route quan- 
tum information in making a modular clocked architec- 
ture. The goal of the router construction is to design a 
system such that application of a field across different 
portions of the device can be used to steer the quantum 
information conditional on where the field has been ap- 
plied. The basic method for achieving this is the gadget 
shown in Fig. 10. The gadget shown can be used to pro- 
duce a state, conditional on where in the device a field 
is applied. Initially the Hamiltonian for this gadget is a 
sum over graph stabilizers for all vertices: 

Hi = -(Z 2 X 1 - Z 2 X V + Z 2 X + Z 1 Z V Z X 2 ) (25) 

Then this Hamiltonian is adiabatically turned off while 
either of the two following Hamiltonians is turned on: 

H f = -(X 1 +X 2 ) (26) 

or 

H' f = -(X v + X 2 ) . (27) 

The point of this is that depending on which of these is 
turned on, the state of the qubit o will differ between 
being a +1 eigenstate of X and a —1 eigenstate of X. 

To see this, first consider the case of ending in Hf. 
Define a stabilizer code using the following operators 

Si = X\Xy , S 2 = X\X , 

Z\ = Z\Z\iZ Q X 2 , X\ = Xi, 

Z 2 = Z 2 Xi : X 2 = X 2 (28) 

where Si and S 2 are the stabilizers for the code and Pj 
are the encoded Pauli operators for this code. Then we 
see that we can express the initial Hamiltonian as 

H l = -(Z 1 +Z 2 (I-S 1 +S 2 )) (29) 

and the final Hamiltonian as 

H f = -(X 1 +X 2 ) (30) 
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FIG. 10. Gadgets and Router. Here we show the stages in 
understanding a router. In (A) we show the basic gadget that 
can be used to conditionally prepare a ±1 eigenstate of X Q 
depending on where the field is applied. (B) is the same con- 
struction, but now demonstrating the flexibility of the con- 
struction to the applied field. The applied field need only 
cover the desired leg in order to get the correct conditional 
behavior. (C) is a non-optimized router construction. If we 
apply the field everywhere but one of the two circled regions, 
this will either route quantum information incoming at 1 and 
2 to 1' and 2' (identity) or 2' and 1' (swap) respectively. Here 
the black nodes are twisted by ^, the grey nodes by — ^ and 
the gradient-shaded node by ty. 

Since S\ and 62 commute with these Hamiltonians, the 
subspaces defined by the eigenvalues of these operators 
will remain constant. The initial Hamiltonian will have a 
ground state which is in the —1 eigenstate of Si and the 
+1 eigestate of S2. Initially the system will also be in the 
+1 eigenvalue eigenspace of Z\ and Z 2 . Upon applica- 
tion of the fields as represented by Hf, the information 
in these later two encoded qubits will be adiabatically 
dragged (with no energy level crossings) to +1 eigenval- 
ues of Xi and X 2 . Thus at the end of this evolution the 
system will be in the +1 eigenvalue eigenspaces of S2, 
X\ and X 2l while being in the —1 eigenvalue eigenspace 
of Si. In particular, X Q — X\S 2 , which implies that the 
system is in the eigenstate of X Q with eigenvalue +1. 

Contrariwise, consider the case ending in Hj. Define 
the code similar to above 

S[ = XiXy,S' 2 = XyX Q , 

Z[ = Z\Z\i Z Q X 2 ,X[ — X\i , 
Z 2 = Z 2 Xy , X 2 = X 2 . (31) 
Then we can express the initial Hamiltonian as 

H i = -(Zi+Z' 2 (-I + S' 1 +S' 2 )) (32) 



and the final, primed, Hamiltonian as 

H' f = -(X[+X' 2 ). (33) 

From these expressions one can deduce that the system 
will start in the —1 eigenstate of the Z 2l S(, and S2 op- 
erators and in the +1 eigenstate of Z' x . At the end of the 
evolution the system will end up in the +1 eigenstate of 
X[ and X' 2 and in the —1 eigenstate of S(, and S' 2 . Since 
X Q = X[S f 2 , this implies that at the end of this adiabatic 
evolution the system will be in the —1 eigenstate of X Q . 

Thus we see that depending on whether the final 
Hamiltonian is Hf or Hp the qubit located at o is in 
either the +1 or —1 eigenstate of X Q . Note that this de- 
pended only on where the field Hf was applied. Further 
note that there is flexibility in spatially achieving this. 
For example, in Fig. 10(B) we show a larger version of 
this gadget. Depending on whether the applied field is 
in either of the two circles, this produces a +1 or — 1 
eigenstate of the last qubit at o. The exact location of 
this applied field is not important except for the point 
that the field spans an entire one of the two "legs" in the 
construction. 

One can extend this idea to then create a router where 
the routing depends only upon where a field has been ap- 
plied. The basic idea is rather simple: if one has the abil- 
ity to conditionally create one of two orthogonal states 
by the location that the field is applied, then one can use 
this as input into a controlled swap gate. In Fig. 10(C) 
we show, for example, such a construction (no attempt 
has been made to optimize this construction) based upon 
the conditional swap gate construction of Smolin and Di- 
Vincenzo [44]. If one applies a field everywhere but one 
of the two circled regions, then this routes the quantum 
information depending on which of the two circled re- 
gions is left out. In particular, information incoming to 
1 and 2 is thus routed (permuted) to output information 
at V and 2'. 



VI. CONCLUSION 

In conclusion, we have introduced a new method for 
building a quantum computer based upon the notion of 
an adiabatic quantum transistor. Two notable benefits 
of this method are that the system is robust to timing er- 
rors (as in universal adiabatic QC) and that it is modular 
in nature (something which prior universal adiabatic QC 
models lacked.) Instead of requiring increasingly accu- 
rate timing and control mechanisms, this model requires 
one to focus on increasing the fabrication quality of en- 
gineered interactions in many-qubit systems. We have 
argued that the noise model for our scheme will follow 
an independent error model and thus is amenable to sta- 
bilization by standard methods of fault-tolerant quantum 
computing. While our actual construct is not optimized 
for current experimental implementation, the mere exis- 
tence of devices like the one we describe, combined with 
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recent experimental progress in building highly control- 
lable quantum simulators [45], gives us hope that adia- 
batic quantum transistors are a viable new path toward 
building a large-scale quantum computer. 

ACKNOWLEDGMENTS 

D.B. was supported by the NSF under grants 0803478, 
0829937, and 0916400. D.B. and G.M.C were sup- 



ported by DARPA under QuEST grant FA-9550-09-1- 
0044. S.T.F. was supported by the IARPA MQCO pro- 
gram. G.M.C. was supported by a Computational Sci- 
ence Graduate Fellowship under DoE grant DE-FG02- 
97ER25308. D.B. thanks Nature and Science referees for 
their career advice. 



[1] W. Shockley, Circuit Element Utilizing Semiconductor 
Material, U.S. Patent 2569347 (Bell Telephone Labora- 
tories, 1948). 

[2] P. W. Shor, in Proceedings of the 35th Annual Sympo- 
sium on the Foundations of Computer Science, edited by 
S. Goldwasser (IEEE Computer Society, Los Alamitos, 
CA, 1994) pp. 124-134. 

[3] D. Aharonov and M. Ben-Or, in Proceedings of the 
29th Annual ACM Symposium on Theory of Computing 
(ACM Press, 1997) pp. 176-188. 

[4] E. Knill, R. Laflamme, and W. H. Zurek, Science 279, 
342 (1998). 

[5] E. Knill, R. Laflamme, and W. H. Zurek, Proc. Roy. Soc. 

London Ser. A 454, 365 (1998). 
[6] D. Deutsch, Proc. Roy. Soc. London Ser. A 425, 73 

(1989). 

[7] S. Bose, Phys. Rev. Lett. 91, 207901 (2003). 
[8] E. Knill, R. Laflamme, and G. J. Milburn, Nature 409, 
46 (2001). 

[9] R. Raussendorf and H. J. Briegel, Phys. Rev. Lett. 86, 
5188 (2001). 

[10] D. Aharonov, W. van Dam, J. Kempe, Z. Landau, 
S. Lloyd, and O. Regev, in 45th Annual IEEE Sympo- 
sium on Foundations of Computer Science (IEEE Com- 
puter Society, Los Alamitos, CA, USA, 2004) pp. 42-51. 

[11] J. Kempe, A. Kitaev, and O. Regev, SIAM Journal of 
Computing 35, 1070 (2006). 

[12] A. Mizel, D. A. Lidar, and M. W. Mitchell, Phys. Rev. 
Lett. 99, 070502 (2007). 

[13] D. Kult, J. Aberg, and E. Sjoqvist, Phys. Rev. A 74, 
022106 (2006). 

[14] W. Wootters and W. Zurek, Nature 299, 802 (1982). 

[15] D. Bacon and S. T. Flammia, Phys. Rev. Lett. 103, 
120504 (2009). 

[16] D. Bacon and S. T. Flammia, Phys. Rev. A 82, 

030303(R) (2010). 
[17] O. Oreshkov, T. A. Brun, and D. A. Lidar, Phys. Rev. 

Lett. 102, 070502 (2009). 
[18] O. Oreshkov, Phys. Rev. Lett. 103, 090502 (2009). 
[19] H. J. Briegel and R. Raussendorf, Phys. Rev. Lett. 86, 

910 (2001). 

[20] S. D. Bartlett and T. Rudolph, Phys. Rev. A 74, 040302 
(2006). 

[21] R. Oliveira and B. Terhal, Quantum Information & Com- 
putation 8, 0900 (2008). 

[22] G. Schaller, S. Mostame, and R. Schiitzhold, Phys. Rev. 
A 73, 062307 (2006). 



[23] A. Doherty and S. Bartlett, Phys. Rev. Lett. 103, 020506 
(2009). 

[24] J. K. Pachos and M. B. Plenio, Phys. Rev. Lett. 93, 
056402 (2004). 

[25] E. Lieb, T. Schultz, and D. Mattis, Annals Phys. 16, 
407 (1961). 

[26] F. Verstraete, D. Porras, and J. I. Cirac, Phys. Rev. 

Lett. 93, 227205 (2004). 
[27] S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 
[28] U. Schollwock, Rev. Mod. Phys. 77, 259 (2005). 
[29] J. Dorier, F. Becca, and F. Mila, Phys. Rev. B 72, 

024448 (2005). 

[30] G. A. Paz-Silva, G. K. Brennen, and J. Twamley, Phys. 

Rev. Lett. 105, 100501 (2010). 
[31] A. Childs, E. Farhi, and J. Preskill, Phys. Rev. A 65, 

012322 (2001). 
[32] P. Pfeuty, Annals of Physics 57, 79 (1970). 
[33] J. Roland and N. J. Cerf, Phys. Rev. A 65, 042308 (2002). 
[34] W. Son, L. Amico, and V. Vedral, Quant. Inform. Proc. 

, 1 (2011). 

[35] N. Schuch, D. Perez-Garcfa, and I. Cirac, Phys. Rev. B 

84, 165139 (2011). 
[36] X. Chen, Z.-C. Gu, Z.-X. Liu, and X.-G. Wen, (2011), 

arXiv:1106.4772 [cond-mat.str-el]. 
[37] S. Bravyi, M. Hastings, and S. Michalakis, J. Math. 

Phys. 51, 093512 (2010), arXiv:: 1001. 0344 [quant-ph]. 
[38] S. Michalakis and J. Pytel, (2011), arXiv:1109.1588 

[quant-ph] . 

[39] D. V. Else, I. Schwarz, S. D. Bartlett, and A. C. Doherty, 

Phys. Rev. Lett. 108, 240505 (2012). 
[40] A. Miyake, Phys. Rev. Lett. 105, 040501 (2010). 
[41] G. Brennen and A. Miyake, Phys. Rev. Lett. 101, 010502 

(2008). 

[42] S. D. Bartlett, G. K. Brennen, A. Miyake, and J. Renes, 

Phys. Rev. Lett. 105, 110502 (2010). 
[43] J. M. Renes, A. Miyake, G. K. Brennen, and S. D. 

Bartlett, (2011), arXiv:: 1103.5076 [quant-ph]. 
[44] J. A. Smolin and D. P. DiVincenzo, Phys. Rev. A 53, 

2855 (1996). 

[45] I. Buluta and F. Nori, Science 326, 108 (2009). 
[46] P Jordan and E. Wigner, Z. Physik 47, 631 (1928). 



16 



Appendix A: Transverse Ising Model with Boundary 
Term Spectrum 

The relevant model on I qubits is 
H,(8) = -(1 - s) Y^[Z]i -sflXl + f^[XUX] i+1 ) . 



i=l 



(Al) 



Because of the single qubit term on the first qubit, both 
the initial and final Hamiltonian are non-degenerate. It 
is convenient to add an extra qubit, call it the Oth qubit, 
and consider the Hamiltonian 



H[{s) = -(1 - a) Y}Z\i -*( E^'^Hi (A2) 



l-l 



i=0 



Then, when we restrict to the +1 eigenspace of [X]q 
([X]o commutes with the Hamiltonian) we will obtain 
Hi(s). Also note that if you conjugate Hi(s) by rii=i[^]i 
you will obtain Hi(s) with the [X]\ term flipped in 
sign. Therefore the —1 eigenvalue of [X]q eigenspace has 
the exact same spectrum as the +1 eigenvalue of [X]q 
eigenspace. Thus we know that H^s) will have exactly 
the same spectrum as Hi(s) but will be two-fold degen- 
erate, the degeneracy corresponding to the eigenvalue of 
[X]q. Thus, if we are interested in the gap of Hi(s) we 
can equally well work with H' t (s) which we now assume. 

The computation of the energy spectrum of H f t (s) fol- 
lows the techniques of Lieb, Schultz and Mattis [25]. Af- 
ter a Jordan- Wigner transform [46] the Hamiltonian can 
be written as 



l-i 



H' l (s)=-sJ2(4-c i )(4 +1 +c i+1 ) 
i 

-(1- S )^(2cjc 2 -/) (A3) 



i=l 



We can express this as 
i 



H[{s) = ^ [Ai,jc\ c j + \{BiA c ) + Bit** 



i,j=0 



TI 



(A4) 



with A a (Z + 1) x (I + 1) symmetric matrix given by 

Aij = -s(S jii+1 + Sij+x) - 2(1 - s)Sij(l - £ i>0 ), (A5) 

B a (I + 1) x (I + 1) anti-symmetric matrix given by 

Bij = -s(5 j:i+1 - Sij+i), (A6) 

and r = (1 — s)l. Then we can find new fermion opera- 
tors, r]k which are linear combinations of the c i and c\s 
such that 

/ / x i x i \ 

H[{s) = ^Ukvlv k + U E A* ~ 2^ Uk + r T 

k=0 V i=0 k=0 ) 

(A7) 



or after simplifying, 

1 ( 1 

H 'i( s ) = (vtVk ~ 2 

k=0 v 



(A8) 



with the cj/cS being the square roots of the eigenvalues of 
4M [25], where 



M 



1 



(A + B)(A-B) 



(A9) 



Note the second Hamiltonian follows by direct calcula- 
tion, but also due to the traceless nature of H f t (s). Ex- 
plicitly, we find that 



Mi, 



s 2 S ZjJ (l - 5 hl ) + (1 - s) 2 5 hJ (l - 5 Zj0 ) 
+ s(l - s)(5 iJ+1 + %+i) . (A10) 



It is easy to check that the (unnormalized) vector 



!*> = £ 



i=0 



l-l 



(All) 



is an eigenvector of M with eigenvalue for all values 
of s. This implies that there is always a mode with zero 
energy; this is exactly a consequence of the fact that the 
Hamiltonian is always two- fold degenerate. As argued 
above we can ignore this fact. 

We can find the eigenvalues of M by the use of an 
ansatz of the form 



7) = 5^sin[fl(i+|)+0]|i). (A12) 



If we apply M to this vector then we obtain three equa- 
tions: two boundary terms and a term from the bulk. 
Define 7 = s(l — s). The bulk term gives us the equation 



7 sin 



(i - |) + + (1 - 2 7 - A) sin 0(i + \) + 



7 sin 



K i + I) + <P 



0, 



(A13) 



where A is the eigenvalue. Some tedious math later, this 
can be turned into 



A = (1 -27) + 27 cos 0. 
The two boundary terms give the equations 

s 2 sin(f +0) +7sin(f +(/))= Asin(f +< 
and 



(A14) 
(A15) 



7 sin 



+ (l-,s) 2 sin 0(1 



A sin 



This first equation can be rearranged to yield 
(1 - s) sin(f +</>)= ssin(f - <j>) 



(A16) 



(A17) 
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while the second one can be manipulated to become 



(l-s)sin 0(i+§)+, 



(A18) 



Solving these equations for ^—^ we obtain the equation 



sm 



+ + 



sin(f+0) sin 6>(Z+ §) +(/> 
This can be reduced to 

cos[6»(Z + 2)] = cos(0Z), 



which has solutions for 



irk 

1 + 1 



(A19) 



(A20) 



(A21) 



where k is an integer. When k = 0, this equation does 
not have a solution for general s. Fortunately, we already 
have the k = eigenvector, so we have found the rele- 
vant non-zero eigenvalues. For k ^ 0, the phase shift (j) 
can be found by solving the transcendental equation of 
Eq. (A17). 

Thus we see that the eigenvalues are 



(s) 2 = 4(l - 2s(l - s) 1 - cos(^) 



(A22) 



for k = 1, 2, . . . , I. Each of these is minimized for s = 1/2, 
for which the eigenvalues become 



(l/2) 2 = 2 + 2cos(^) =4cos 2 (^). (A23) 



From this we see that u^s smallest value is 0(1//), oc- 
curring when k = I. 



